Necessary Optimality Conditions for Fractional Action-Like 
Problems with Intrinsic and Observer Times 
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Abstract: We prove higher-order Euler-Lagrange and DuBois-Reymond stationary conditions to fractional action- 
like variational problems. More general fractional action-like optimal control problems are also considered. 
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1 Introduction 

The study of fractional problems of the calculus of 
variations and respective Euler-Lagrange type equa- 
tions is a subject of strong current research because 
of its numerous applications: see e.g. HI [3l 01 
H US H21 El HI- F Riewe ID] obtained a ver- 
sion of the Euler-Lagrange equations for problems of 
the calculus of variations with fractional derivatives, 
that combines both conservative and non-conservative 
cases. In 2002 O. P. Agrawal proved a formula- 
tion for variational problems with right and left frac- 
tional derivatives in the Riemann-Liouville sense Ofl. 
Then, these Euler-Lagrange equations were used by 
D. Baleanu and T. Avkar to investigate problems with 
Lagrangians which are linear on the velocities Q. 
In |[T2l fractional problems of the calculus of varia- 
tions with symmetric fractional derivatives are con- 
sidered and correspondent Euler-Lagrange equations 
obtained, using both Lagrangian and Hamiltonian for- 
malisms. In all the above mentioned studies, Euler- 
Lagrange equations depend on left and right fractional 
derivatives, even when the problem depend only on 
one type of them. In lfl3ll problems depending on 
symmetric derivatives are considered for which Euler- 
Lagrange equations include only the derivatives that 
appear in the formulation of the problem. In (31 [H 
Euler-Lagrange type equations for problems of the 
calculus of variations which depend on the Riemann- 
Liouville derivatives of order (a, (3), a > 0, j3 > 0, 
are studied. 

In (51 HH [I7j, C. Udriste and his coauthors re- 
mark that the standard multi-variable variational cal- 
culus has some limitations which the multi-time con- 
trol theory successfully overcomes. For instance, the 



classical multi-variable variational calculus cannot be 
applied directly to create a multi-time maximum prin- 
ciple. In (6l 13 two-time Riemann-Liouville frac- 
tional integral functionals, depending on a parame- 
ter a but not on fractional-order derivatives of or- 
der a, are introduced and respective fractional Euler- 
Lagrange type equations obtained. In [03, Jumarie 
uses the variational calculus of fractional order to de- 
rive an Hamilton- Jacobi equation, and a Lagrangian 
variational approach to the optimal control of one- 
dimensional fractional dynamics with fractional cost 
function. Here, we extend the Euler-Lagrange equa- 
tions of J6l [7J by considering two-time fractional 
action-like variational problems with higher-order 
derivatives. A DuBois-Reymond stationary condition 
is also proved for such problems. Finally, we study 
more general two-time optimal control type problems. 

2 Preliminaries 

In 2005, El-Nabulsi (cf. 0) introduced the following 
Fractional Action-Like VAriational (FALVA) prob- 
lem. 

Problem 1 Find the stationary values of the integral 
functional 



/[?(•)] 



1 



r(a) J a 



L(e, q (6),q(e))(t-er- 1 de 



(Pi) 



under given initial conditions q(a) = q a , where q = 
■S, T is the Euler gamma function, < a < 1, 9 is 
the intrinsic time, t is the observer time, t ^ 0, and the 



Lagrangian L : [a, b] x M. n x M. n - 
with respect to all its arguments. 



. is a C 2 -function 



Along all the work, we denote by d{L the partial 
derivative of a function L with respect to its 2-th argu- 
ment, ieN. 

Next theorem summarizes the main result of |[6l . 

Theorem 2 (cf. f^jl)Ifq(-) is a solution of Problem^ 
that is, q(-) offers a stationary value to functional 
<\Pity, then q(-) satisfies the following Euler-Lagrange 



equations: 



d 2 L(9,q(9),q(9)) - -d 3 L(9,q(9),q(t)) 



ft 



d 3 L(e,q(6),q(9)) . (1) 



In this work we begin by generalizing the Euler- 
Lagrange equations (Q]) for FALVA problems with 
higher-order derivatives. 

3 Main Results 

In $3J] and we study FALVA problems with 
higher-order derivatives. The results are: Euler- 
Lagrange equations (Theorem [5]) and a DuBois- 
Reymond stationary condition (Theorem ITQb for such 
problems. Then, on section £|3.3[ we introduce the 
two-time optimal control FALVA problem, obtaining 
more general stationary conditions (Theorems [TBI and 

\m. 

3.1 Euler-Lagrange equations for higher- 
order FALVA problems 

We prove Euler-Lagrange equations to higher-order 
problems of the calculus of variations with fractional 
integrals of Riemann-Liouville, i.e. to FALVA prob- 
lems with higher-order derivatives. 

Problem 3 The higher-order FALVA problem con- 
sists to find stationary values of an integral functional 



1 



L[9,q(9),q(9), 



r(a) Ja 

...,q^\e)){t-er- 1 de, (p m ) 

m > 1, subject to initial conditions 



q (t) (a) 



0, . . . , m — 1 , 



(2) 



where q(°\0) = q(9), q <y% \9) is the i-th derivative, 
i > 1; r is the Euler gamma function; < a < 1; 9 is 
the intrinsic time; t the observer's time, t ^ 9; and the 
Lagrangian L : [a, b] x IR nx ( m+1 ) — > R is a function 
of class C 2m with respect to all the arguments. 



Remark 4 In the particular case when m = 1, func- 
tional \P m ) reduces to \P\\ and Problem\3\to\l\ 



To establish the Euler-Lagrange stationary condi- 



tion for Problem flP m ] ), we follow the standard steps 
used to derive the necessary conditions in the calculus 
of variations. 

Let us suppose g(-) a solution to Problem [3] The 
variation 5I m [q(-)] of the integral functional \P m ) is 
given by 

where 8q^> € C 2m ([a, b\; M. n ) represents the varia- 
tion of q( l \ i = 1, . . . , to, and satisfy 

5q (i) (a)=0. (4) 

Having in account conditions (O, repeated inte- 
gration by parts of each integral containing 5q^ in 
© leads to 



to = 1 : SI[q(-)} = 

1 — ft 
~ t-9 

m = 2: 5I 2 [q(-)]- 



r(a) J a 



d 2 L - ±D>L 



(t - 9) a ~ l ■ 5qd9 ; (5) 



1 



r(a) Ja 



d9 



;i -ft)(2-a) 



d 4 L 



{t-9) 2 1 
and, in general, 

SI m [q(-)] 




r(ft 
i — a 



{t - 9) a ~ l ■ 5qd9 ; (6) 



t 



i=l 



in k 



srsr , r(»-a + i) / k \ a*-* 



• (t - 9) a ~ l ■ 5qd9 . 

The integral functional I m [-) has, by hypothesis, 
a stationary value for q(-), so that 

5I m [q(-)]=0. 

The fundamental lemma of the calculus of variations 
asserts that all the coefficients of 5q must vanish. 



Theorem 5 (higher-order Euler-Lagrange equa- 
tions) If q(-) gives a stationary value to func- 
tional < \P m h then g(-) satisfy the higher-order Euler- 
Lagrange equations 



EM)' 



i=0 



di+2L 



l (9)A(9),..^q {m) {9) 



where m > 1 and 

F(e, q (e)A{o),...,q( 2m - x Xe)) 



(V) 



Proof: Theorem[5]is proved by induction. For m = 1 
and m = 2, the Euler-Lagrange equations CD and (©- 
(TTOb are obtained applying the fundamental lemma of 
the calculus of variations respectively to (f5]) and ©. 
From the induction hypothesis, 



i=0 



1 — a 



ji—l 



j=l 



+ED-D 



.! r(i-a + l) 



fc=2 8=2 



7 9k+2L 



(t - eyr(i -a)\k-i) d9 k - i k+2 



with the partial derivatives of the Lagrangian L eval- 
uated at (9,q{e),q(9),...,q( m \e)). 

Remark 6 Function F in ^} may be viewed as an 
external non-conservative friction force acting on the 
system. If a = 1, then F = and equation (0 is noth- 
ing more than the standard Euler-Lagrange equation 
for the classical problem of the calculus of variations 
with higher-order derivatives: 

£(-!)'— $ +2 L (e, q(9),q(9), - , 1 (m) (*)) = . 
i=0 



Remark 7 If m = 1, the Euler-Lagrange equations 
((TJ) coincide with the Euler-Lagrange equations (Q]). 

Remark 8 For m = 2, the Euler-Lagrange equations 
© reduce to 

d 2 L (9, q, q, q) - ^<9 3 L (9, q, q, q) 



d 2 

+ 
where 

F(9,q,q,q,q 



d Q2 d ^ L ^ 



1 — a 



F(9,q,q,q,q) (9) 



d 3 L(9,q,q,< 



2—d 4 L(8,q,q, 



r(3 - a) 



(i - #) 2 r(l - a) VO 



1 — a 



d 3 L (9, q, q, q) - 
(1 -q)(2-a) 



d±L{9,q,q,q) 
2 -^d 4 L (9, q, q, q) 



d 4 L(9,q,q,q) . (10) 



(8) 



= F (fl, g(0), . . . , gW- 1 )^)) , m = j > 2. 

(11) 

We need to prove that equations ©-([U) hold for m = 
j + 1. For simplicity, let us focus our attention on the 
variation of q^ +l \9). From hypotheses (variation of 
q^ +1 \9) up to order m = j), and having in mind 
that C{ + Cf +1 = C(+l and mT{m) = T(m + 1), we 
obtain equations © for m = j+l using integration by 
parts followed by the application of the fundamental 
lemma of the calculus of variations. □ 

It is convenient to introduce the following quan- 
tity (cf. 03]): 



m-j 



d 



& = E (- 1 )^^+i+2^ (9, q(9),q(9), q^ (9) 



i=0 



(12) 

j = 1, . . . , m. This notation is useful for our purposes 
because of the following property: 

^ = d j+1 L (9, q(9), q(9), q^ (9)) - 

(13) 

j = l,...,m. 



Remark 9 Equation (0 can fee written in the follow- 
ing form: 



n,L[9,q(9),q(9),---,Q {m) (0))-^ 1 



F 



(9,q(9),q(9),...,q( 2m -V(8)) . (14) 



3.2 DuBois-Reymond condition for higher- 
order FALVA problems 

We now prove a DuBois-Reymond condition for 
FALVA problems. The result seems to be new even 
for m = 1 (Corollary ITU). 

Theorem 10 (higher-order DuBois-Reymond con- 
dition) A necessary condition for g(-) to be a solution 
to Problem \3\ is given by the following higher-order 



DuBois-Reymond condition: 



-{L[9,q(9),q(9),- <1 



(rn 



3=1 



= d 1 L(e,q(e),q(e),..., q ( m \9) / 

+ F(9, q (6), q (6),...,q( 2m - 1 \e))- q (6), (15) 

where F and ip 1 are defined by ([8]) and (|12l) . respec- 
tively. 

Remark W If a = 1, then F = a?i<i condition 
(1151 ) is reduced to the classical higher-order DuBois- 
Reymond condition (see e.g. M5\l ): 



d\L (o, q(6),q(9), . . . , q^ m 



^ { L[0.q[9).,ii9) u 



(m 



Proof: The total derivative of 

L(9,q(9),q(9),---,<l {m 
with respect to 9 is: 
d_ 

1e 



3=1 



^■q^(9) 

3=1 



L(9,q(9),q(9),...,q 



(m 



E^'^V) 



3=1 



^(9,q(9),q(9),...,q^(9)) 



m 

- E ' * U) ( e ) + ^ ■ • (16) 

3=1 

From (fT3T ) it follows that (fT6l ) is equivalent to 



d0 



L(e,q{9),m,...,q^ m 



3=1 



dL 

~d9 



i H \^ f i H \^ ) i * * * 1 



Jrri 



dL 



dqti) 

m 

Elv^iy-U 



^(0),g(0),...,g( m )(0)) 



(17) 



We now simplify the last term on the right-hand side 
of <T7): 



E 



m— 1 

E[ 

j'=o 



<9L 



Z (j+ 1 )(0)-^'.g(i+ 1 )(0) 



+ ^J+1 . g O"+2)(0) 



m— 1 

E 

i=o 



<9L 



,C?+1)(0) 



-v; -g(fl)+v m v m+1) (#) 



(18) 



where the partial derivatives of the Lagrangian L are 
evaluated at (9,q(9),q(9), . . . , q^ (9)). Substitut- 



ing (1181) into (U7I ). and using the higher-order Euler- 
Lagrange equations (0, we obtain the intended result, 
that is, 



d_ 

7w 



L 



...,q 



(m 



E^'-^w 



^(9,q(9),q(9),...,q( m \9)) 

}, q (9),q(9),...,q( m \9)) -^(9) 
+ V° -g(0) -4> m ■q {m+1) {9) 

d 1 L(9,q(9),q(9),...,q( m 



+ 



dL 

dq( m ) 



+ F 



•> "i V / J ^ V /' ' ' ' ' 



(2m- 



since, by definition, 
<9L 



0gM 



and 



V>° = E^ 1 )^^ (fl, ^(0), g((9), ...,<? 



i=0 



(m 



□ 



Corollary 12 (DuBois-Reymond condition) Ifq(-) is 
a solution of Problem [7] then the following (first- 
order) DuBois-Reymond condition holds: 

^ {L (9, q(9),q(9)) - d 3 L (9, q(9),q(9)) ■ q(9)} 
= d\L (9, q (e), q (9))+\^d 3 L (9, q(9),q(9))-q(9) . 



(19) 



Proof: For m = 1, condition (fT5l) is reduced to 



is equivalent to 



(16 



{L (0, q(8),q(8)) — ip 1 ■ q(0)] 

d\L (8, q(6), q(8)) + F (8, q(8),q(8)) ■ q(8) . 

(20) 



1 



Having in mind ([8]) and (fl2l) . we obtain that 

^ = d 3 L(8,q(8),q(8)) , 



(21) 



F(8,q(8),q(8)) 



1 — a 



d 3 L(8,q(8),q(8)) . (22) 



One finds the intended equality (fl9l by substituting 
the quantities (EJ and (|22]> into ([20]). □ 



3.3 Stationary conditions for optimal control 
FALVA problems 

Fractional optimal control problems have been stud- 
ied in (21 [H|9l. Here we obtain stationary conditions 
for two-time FALVA problems of optimal control. We 
begin by defining the problem. 

Problem 13 The two-time optimal control FALVA 
problem consists in finding the stationary values of the 
integral functional 



i 



r(a) J a 

when subject to the control system 



L(8,q(8),u(8)){t-8) a - 1 d9 1 
(23) 



q(9)=ip(8,q(8),u(8)) 



(24) 



and the initial condition q(a) = q a . The Lagrangian 
L : [a, t] x 1" x W —> R and the velocity vector 
(p : [a, b] x R n x f -> 1" are assumed to be C 1 
functions with respect to all their arguments. In ac- 
cordance with the calculus of variations, we suppose 
that the control functions u(-) take values on an open 
set ofW. 

Remark 14 Problem \J} is a particular case of Prob- 
lem |72] where (p(6, q, u) = u. FALVA problems of 
the calculus of variations with higher-order deriva- 
tives are also easily written in the optimal control 
form (123l )- (|24l) . For example, the integral functional 
of the second-order FALVA problem of the calculus of 
variations, 



1 



r(a) J a 



L(8,q(8),q(6),Q(0))(t-Or- l d8, 



r(a) J a 



L(8,q°(8),q\8), U (8))(t-8) a - 1 d8, 

U{8)=q\8), 
\q\8) = u{8). 



We now adopt the Hamiltonian formalism. We 
reduce (T23l-(l24]) to the form \P\\ by considering the 
augmented functional: 

1 rt 



r(a) Ja 



[H(8,q(8),u(8),p(8)) - p(8) ■ q(8)]d8 , 



(25) 



where the Hamiltonian H is defined by 

n(8,q,u,p)=L(8,q,u) {t-8)^ 1 +p-cp(0,q,u) . 

(26) 

Definition 15 (Process) A pair (<?(•), u(-)) that satis- 
fies the control system q(8) = <p (8 , q(8) , u(8)) and 
the initial condition q[a) = q a of Problem\13\is said 
to be a process. 

Next theorem gives the weak Pontryagin maxi- 
mum principle for Problem [T3l 

Theorem 16 If (q(-),u(-)) is a stationary process for 
Problem \T3\ then there exists a vectorial function 
p(-) G C 1 ([a,6];R n ) such that for all 8 the tuple 
(q(-),u(-),p(-)) satisfy the following conditions: 

• the Hamiltonian system 

(q(8)=d 4 H(8,q(8),u(8),p(9)), 



p(8) = -d 2 H(8,q(8),u(8),p(8)) 



(27) 



the stationary condition 

d 3 H(8,q(8),u(8),p(8))=0; (28) 



where Ti is given by (126I ). 

Proof: We begin by remarking that the first equation 
in the Hamiltonian system, q = d^H, is nothing more 
than the control system (|24l . We write the augmented 
functional (T2B1 in the following form: 



1 



r(«) Ja 



H-p(8)-q(8) 
(t - 8)*- 1 



(t - 8) a ~'d8 , (29) 



where TL is evaluated at (9, q(9),u(9),p(9)). Intended 
conditions are obtained by applying the stationary 
condition CD to ( f29l ): 



( d d 


H-pq 


d 


H-p-q 1 


i-a a 


H-pq 


| de dq 


(t-e)'*- 1 


~ dq 




t-e dq 


(t-e)"- 1 


1 d a 


H-p-q 


_ d_ 


H-pq 


1-as a 


T w-p-9 1 


^ de dii 


(t-e)"- 1 


du 


(t-e) *- 1 


t-0 9m 





-p" = <9 2 w 
o = s 3 h 



□ 

Remark 17 For problems of the calculus of 

variations, Theorem 1 7 61 tafces the form of Theorem\5\ 

Definition 18 (Pontryagin FALVA extremal) We call 
any tuple (q(-),u(-),p(-)) satisfying Theorem \16\ a 
Pontryagin FALVA extremal. 

Next theorem generalizes the DuBois-Reymond 
condition dT3T > to Problem [TBI 

Theorem 19 The following property holds along the 
Pontryagin FALVA extremals: 

^(9, q(9),u(9),p(9)) = 8^(9, q(9),u(9),p(9)) . 



Proof: Equality ( f30T > is a simple consequence of The- 
orem [16] □ 

Remark 20 In the classical framework, i.e. for a = 
1, the Hamiltonian TL does not depend explicitly on 9 
when the Lagrangian L and the velocity vector tp are 
autonomous. In that case, it follows from (1301 ) that the 
Hamiltonian TL (interpreted as energy in mechanics) 
is conserved. In the FALVA setting, i.e. for 1, this 
is no longer true: equality (1301) holds hut we have no 
conservation of energy since, by definition (cf. (1261 )), 
the Hamiltonian TL is never autonomous (TL always 
depend explicitly on 9 for a/1, thus d\TL ^ 0). 
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